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Abstract 



A scenario for systems with slow dynamics is characterised by stating 
that there are several temperatures coexisting in the sample, with a single 
temperature shared by all observables at each (widely separate) time-scale. 

In preparation for the study of granular rheology, we show within this 
framework that glassy systems with driving and friction that are generic and 
do not correspond to a thermal bath — and whose microscopic 'fast' motion is 
hence not thermal — have a well-defined macroscopic temperature associated 
to the slow degrees of freedom. 

This temperature is what a thermometer coupled to the system will mea- 
sure if tuned to respond to low frequencies, and since it can be related to 
the number of stationary configurations, it is the formalisation of Edwards' 
'compactivity' ideas. 



Granular matter set into motion by shearing, shaking or tapping is one of the most inter- 
esting cases of macroscopic out of equilibrium systems. Given a granular system subjected 
to some form of power input that makes it perform stationary flow on average, a very natural 
question that arises is to what extent it resembles a thermodynamic system of interacting 
particles such as, for example, a liquid. 

More specifically, many attempts have been made to define a 'granular temperature' (see 
e.g. EH). In order to deserve its name, a temperature has to play the role of deciding the 
direction of heat flow: it must be connected to a form of the zero-th law. In order to pursue 
this line, however, one has to somehow take care of the characteristics of granular flow that 
distinguish it from usual kinetic theory: 

1) Energy is not conserved, and, more generally, the motion does not have the very 
strong phase-space volume conservation properties typical of Hamilton's equations. The 
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dissipation is due to friction which is in general not linear in the velocity, and dependent 
upon the relative positions of the particles. 

2) Power is supplied by tapping (which may be periodic in time) or by shearing, a manner 
very different from that of the 'collisions' of a thermal bath. 

Under these circumstances, there is no reason why the observables should be related to 
a Gibbs (or equivalent) ensemble, and the possibility of having thermodynamic concepts 
seems lost. 

In this paper we shall consider situations with shear and friction, in the limit of weak 
shear. The effect of a coherent 'tapping' will be discussed in further work!. In this limit of 
'slow rheology', it will turn out that even though the rapid motion cannot be associated with 
thermal motion, there appears for the slow flow a natural temperature playing the usual role 
in thermometry and thermalisationS. 

The computation presented here can be done in a wide range of approximation schemes 
consisting in resumming a perturbative expression for the dynamics in several forms — and 
to higher and higher levels of approximation — in particular the so-called mode-coupling 
approximation. Here, for concreteness, I will carry it through for a simple model for which 
the (single mode) mode-coupling approximation is exact. 

Multiple Thermalisation in Aging and Rheology 

Granular systems have been recognised as being closely related to glassy systems!!. Ac- 
cordingly, several recent developments and models have been borrowed from the field of 
glasses to understand their properties^. 

A picture has arisen in the last few years for aging or gently driven glasses involving 
multiple thermalisations at widely separated timescales (seeH for a review). In the simplest 
scheme, the situation is as follows: Given any two observables A and B belonging to the 
system, define the correlation function as: 

(A(t)B(t'))=C AB (t,t') (1.1) 

and the response of A to a field conjugate to B: 
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sh B (F) {m = RabM (L2) 

For a pure relaxational (undriven) glass, the correlation breaks up into two parts: 

C AB (t,t') = C F AB {t - t>) + C AB (1.3) 

with h the same growing function for all observables A, B. The fact that CAB{t,t') never 
becomes a function of the time-differences means that the system is forever out of equilib- 
rium, it ages. If instead the glass is gently driven (with driving forces proportional to, say, 
e), aging may stop, and we have: 

C AB (t, = C F AB (t - t') + C AB (—) (1.4) 
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where r G is a time scale that diverges as e goes to zero. 

In the long time limit and in the small drive limit, the time scales become very separate. 
When this happens, it turns out that the responses behave as: 

R AB (t,t')=(3^C^ B + (3*^C AB (1.5) 

in the aging and the driven case. 

The fast degrees of freedom behave as if thermalised at the bath temperature f3. On 
the other hand, the effective, system-dependent temperature T* = 1/(3* indeed deserves its 
name: it can be shownEH] that it is what a 'slow' thermometer measures, and it controls the 
heat flow and the thermalisation of the slow degrees of freedom. It is the same for any two 
observables at a given timescale, whether the system is aging or gently driven. Furthermore, 
it is macroscopic: it remains non-zero in the limit in which the bath temperature is zero. 

If the system is not coupled to a true thermal bath, but energy is supplied by shaking 
and shearing, while it is dissipated by a nonlinear complicated friction, there is no bath 
temperature j3. What will be argued in what follows is that even so, the 'slow' temperature 
(3* survives despite the fact that the fast motion is not thermal in that case. Indeed, if we 
have correlation having fast and slow components: 

C AB (t,t') = C^ B (t,t') + C s AB (t,t') (1.6) 

the response is of the form: 

d_ 

dt' 



R AB (t,t>) = R F AB (t,t')+p*—C s AB (t,t') (1.7) 



with the fast response R AB (t,t') bearing no general relation with the fast correlation 

The effective 'slow' temperature so defined is then found to be directly related to Ed- 
wards' compactivityl, but, in the spirit of Ref.0, in the context of slowly moving rather than 
stationary systems^. It seems also closely related to the macroscopic temperature driving 
activated proceses in the SGR modefi3. 



A Simple Example 

For concreteness, let us consider a variation of the standard mean-field glass modelEli. 
The variables are Xi, i = 1, N, and are suject to an equation of motion: 

mx % + + nx t = - € fl sheaT \x) - /; friction '(i) (1.8) 

dXi 

The left hand side is just Newtonian dynamics (with Q possibly time-dependent), with 
a 'glassy' potential which we can take, for example, as: 

E(x) = JijkXiXjXk (1.9) 

where the is a symmetric tensor of random quenched variables of variance 1/N 2 . These 
terms correspond to the p-spin glass0. It was realised some ten years ago that this kind of 
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model constitutes a mean-field caricature of fragile dass0, and in particular its dynamics 
yields above the glass transition the simplified mode-coupling equations^!. 

On the right hand-side of Eqn. ( |1.8| ) we have added two terms that mimic granular 
experiments. The forces /j do not derive from a potential, for example^!: 

/• Shear; \x) = K%XiX k (1.10) 

with Kj k a non-symmetric tensor with random elements with variance 1/N. They pump 
energy into the system, and hence play a role similar to shearing. All our discussion will 
be restricted to weak driving, i.e. e small. For the friction terms we can take, instead of a 
linear term cx £i, a more complicated odd function y. friction = j friction 

These equations for the correlation C(t,t') = jfJ2{xi{t)xi(t')) and response R(t,t') = 
jj S(xi(t))/Shi(t') can be exactly solved in the large N limit. One may do so by reducing 
the system to a self-consistent single-site equationE-3 

mx + Qx = - e /' shear '(t) - f Mction \x) + 

+ 6 f dt'C(t, t')R{t, t')x(t') + p(t) (1.11) 

/ s ^ ear (t) and pit) are independent coloured Gaussian noises that satisfy: 

< /' shear ' (t)/' shear V) >=< p(tW) >= 3C 2 (t,f) (1.12) 
Equation (|1.11|) is supplemented by the self-consistency conditions 

(x(t)x(t'))=C(t,t') ; ^(M0 = ^|y (1.13) 

where h{t) is a field that acts additively in fll.llQ . 

We now perform the usual step of separating 'fast' and 'slow' functions. Accordingly, we 
put: 

C(t,t')=C F (t,t') + C s (t,t') ; R(t,t') = R F (t,t') + R s (t,t') 
^shear' = ^shear' + ^hear' . p = pp + ps (L14) 

the additive separation is such that in the region of large time-differences Cp{t,t') and the 
integral of Rp(t, t') tend to zero, and the induced noises are now divided into fast and slow: 

< /> hear '(t)/> near '(t') >=< PF it)p F {t') > = 3Cf (MO 

< fs She&T \t)fs Shear \t') >=< Ps(t)ps(t') > = 3Cg(t,0 (1.15) 

We can now rewrite the single-site equation ( |1.11| ) in the following way 

mx + Qx = -e/> near '(t) - /' friction '(x) + 

+ 6 J* dt'C F (t,t')R F (t,t')x(t') + p F (t) + Z(t) (1.16) 
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where: 

Z{t) = -ef s sheaT \t) + 6 fdl/C s (t,lf)R s {t,lf)x(lf)+p s {t) (1.17) 



Consider equation ( |1.16[ ): it describes a single degree of freedom which has nonlinear 



friction and a (short) memory kernel, plus a slowly varying field Z(t). Upon the assumption 
of large separation of timescales (which will be valid if the system is weakly sheared and 
'old'), we can treat Z(t) as adiabatic. However, because of the absence of detailed balance, 
we know that the distribution for fixed Z is non-Gibbsean. The fast correlation and response 
functions will be of the form Cp(t, t') = Cp(t — if) and Rp(t, if) = Rp(t — t'), but we cannot 
say anything in general about their relation. The average of ((x))z over an interval of time 
large compared to the fast relaxation (the range of Cf and Rf) is a certain function of Z: 

«*»* = 8J f£ l (1-18) 

which defines the single variable function F(Z). 

We now turn to the slow evolution. Because the memory kernels in equation (|1.17|) are 



slowly varying and smooth, we can substitute x by its average ((x))z- We hence have: 

Z{t) = -ef s sheai \t) + 6 fdlfC g {t,lf)R 8 (t,e)?^(e)+p s (t) + h a ^ h (t) (1.19) 



where we have explicitated the slow component a field acting additively in ( p.. 11 ). The 
self-consistency equations now read: 

r (t f) - ( 9F{Z) (t) 9F{Z) (t>)\ ■ R (t f) - 6 > dF{Z) (tY. h 20) 

Cs ^ t) - { ^ Z - {t) ^Z- {t)) ' Rs{tjt) - 5h ^diab it ^-dZ-^ (L20) 



Equations (|1.19 ) and ( 1.20 ) describe the slow part of the evolution. The only input of 



the fast equations is through the function F(Z). The manner of solution depends little on 
the fact that the 'fast' evolution is not thermal, and is by now standard. Here I sketch the 
stepsEl for completeness. Our aim is to show that they admit in the small e limit a solution 
of the form: 

where (3* is the effective temperature that emerges for the slow dynamics, to be determined 
by the matching with the 'fast' time-sector. To show this, we first write r = hx(h(t)), 
t' = \n(h(t')), etc, and put 

C S {t,t') = C{T-T') 
Rs(t,t')=f3*-^C(T-T') 

Z(t) = Z(t) (1.22) 
Equations ( |1.19| ) and (|1.20p now take the form, in the e — > limit: 
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Z(t) 



shear' 



(r) + 6/3* / dr'C{r - t')C'{t - t') 



dZ 



(r') + p s (r) + h 



adiab 



1.23) 



where we have explicitated the slow component a field acting additively in ( |1 . 1 1| ) . The 
self-consistency equations read: 



and 



C(r - r') 



ldF{Z) \r) dF{Z \r>)) 



dZ 



dZ 



;i.24) 



5h adiah (T>) dZ 



1.25) 



To prove that ( |1.22| ) is a solution of the system ( |1.23| ) - ( |1.25| ), and hence that ( |1.21| ) is 
a solution of ( |1.19|) and ( |1.20| ) we can proceed as follows: We introduce an infinite set of 
auxiliary variables h(r) and a set of dynamical variables i/i satisfying an ordinary Langevin 
equation with inverse temperature (3*: 



d 2 



d 



9F (gj A m 



with 



Ur)^(r w )) = 2T*T j S lj S(T-r w ) . 



1.26) 



;i.27) 



We can now choose the Aj and the lj such that: 



j 

Y,A J AjH J (r)=h adiah (r) 



(1.28) 



and check that the quantity J2j Aji/j obeys the same equation of motion as Z{r). 

Because the problem reduced to an ordinary (not glassy!) Langevin equation, we can 
assure that the system thermalises at temperature T* = 1/(3*, and hence verify that the 
ansatz closes. 

We have followed essentially the same steps as in the treatment of a glassy models coupled 
to a 'good' heat bath. Here instead of having two time scales each with its own temperature 
we have a temperature associated only with the low frequency motion. 
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Conclusions 

It has been recognised for some time that granular systems bear a deep similarity with 
glassy systems at essentially zero temperature. In order to introduce some agitation, both 
tapping and shearing have been often introduced. This certainly makes sand look more like 
a fluid, but at the same time poses the problem of introducing (and dissipating) energy in 
a manner that is quite different from that of a thermal bath. 

We have shown in this paper that, at least within an approximation scheme and for 
slowly flowing systems, one can still introduce thermodynamic concepts — provided one 
attempts to apply them only to the low frequency motion. 

For very small power input, it turns out that the fluctuation-dissipation temperature of 
the slow degrees of freedom coincides in these models with Edwards', defined on the basis of 
the logarithm of the number of stable configurations!. For stronger driving power, an anal- 
ogous definition needs the counting of stable states each composed of many configurations. 
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